Nonlinear statistical damage constitutive model of granite based on the energy dissipation ratio

The stress–strain curves and mechanical properties of Shuangjiangkou granite were obtained using five groups of conventional triaxial tests under various confining pressures using MTS815 rock test equipment. From the microscale, mesoscale, and macroscale perspectives, four types of mechanisms that contribute to energy dissipation during granite deformation were investigated. Based on the energy dissipation ratio, a new approach for estimating crack closure stress and damage stress is proposed. The energy dissipation ratio was substituted into the Weibull distribution function, and then a new nonlinear statistical damage constitutive model of granite based on the energy dissipation ratio was constructed after Biot’s theory was modified per the Lemaitre strain equivalence principle. By comparing experimental data with theoretical values estimated by the model, the model’s rationality and correctness were confirmed.

The mechanical behaviour and characteristics of rock under diverse external loads have always been a subject of interest to geotechnical engineers [1][2][3][4][5][6][7] . The damage, deformation, failure and post-failure deformation of rock can be directly used for large rock engineering design. The description of rock strength and deformation behaviour (especially the study of constitutive relationships) is the basis for evaluating the stability of rock engineering structures [8][9][10] . In general, the most effective technique to investigate the rock deformation process and failure mechanism is to determine the complete stress-strain relationship. According to a significant number of studies, the complete stress-strain curve of brittle rock can be split into five stages, namely, compression closure of primary cracks, linear elastic deformation, stable crack growth, unstable crack growth, and post-peak failure deformation, by four stress thresholds [11][12][13] . It is critical to determine stress thresholds at each step of rock development to investigate its strength and deformation 14,15 . To date, there have been several approaches for determining stress thresholds, but none for determining the thresholds in terms of the energy evolution law.
The deformation and fracture of rocks is the outcome of an irreversible dissipation process in which energy is converted from one form to another and entropy is generated, according to thermodynamics 16 . Detailed investigation of energy transmission and transformation during the evolution of rock damage can more accurately reflect the law of rock damage. However, in many studies [17][18][19][20] , the mechanism of energy dissipation at the transscale microscopic, mesoscopic, and macroscopic levels of damage and failure progression in brittle rock has not been investigated in great detail. An acoustic emission (AE) is defined as a transient elastic wave generated by the rapid release of energy within a material, therefore, it can be used to predict rock failure as well as to study earthquake processes [21][22][23][24][25][26][27][28] . Brittle rock materials are naturally heterogeneous, with macroscopic discontinuous structures (such as flaws, joints, cracks, and weak surfaces) as well as internal composite microscopic structures (such as different mineral compositions and randomly distributed particles). Even rocks that appear intact contain primitive microcracks and microvoids. Under the action of external factors (including force, temperature, electricity, magnetism, radiation, etc.), a large number of new microfissures and micropores will form inside the rock, and these microfissures will nucleate, cross and merge, resulting in gradual damage and even destruction of the rock 29,30 . In damage mechanics [31][32][33][34][35] , all microdefects are continuous, and their effect on the material is expressed by one or more internal field variables (i.e., damage variable), but these microdefects are discrete and random in nature. www.nature.com/scientificreports/ As a result, the constitutive relationship of rock may be investigated from a statistical perspective, and the constitutive theory and evolution equation of materials with damage can be studied using the random distribution function.
The Weibull distribution is one of the most extensively used statistical models for heterogeneous brittle materials, and the statistical uniformity evaluation approach paired with the Weibull distribution has been widely utilized to analyse brittle material failure [36][37][38][39] . The majority of current studies use the mechanical or deformation properties of rock as the foundation for developing a statistical damage constitutive equation. Krajcinovic et al. 32 used the Weibull distribution function to establish the damage equation for brittle materials and examined the related parameters. Based on the statistical evolution equation, Cao et al. 40 developed a statistical damage constitutive model for rock strain softening and hardening under conventional triaxial compression. Pourhosseini et al. 41 suggested a constitutive model based on the Mohr-Coulomb failure criterion that may account for both rock strain-softening and dilation behaviour. Molladavoodi et al. 42 assumed that the shear and bulk modulus of rock elements follow the Weibull distribution and then used a constitutive model to link the Mohr-Coulomb strain-softening model to a commercial discrete element method (DEM) code. Although these studies used mathematical statistics to develop rock constitutive models, none of them are based on energy evolution rules, which could better describe the nature of rock deformation and fracture. Despite the fact that the rule of energy evolution can in essence describe the process of rock deformation and failure, few studies have included it in the damage constitutive relationship of mesoscopic statistics.
The stress-strain relationship curves for five groups of rock samples were acquired using MTS815 rock triaxial test equipment. This equipment was used to conduct conventional triaxial compression tests of granite under various confining pressures to investigate the relationship between energy dissipation and the damage constitutive equation of Shuangjiangkou granite under different in situ stresses. From the microscale, mesoscale, and macroscale perspectives, four types of energy dissipation mechanisms during rock damage evolution were analysed. A new approach for calculating stress thresholds (crack closure stress and damage stress) based on the energy dissipation ratio is provided. The energy dissipation ratio was substituted into the Weibull distribution function, and generalized Hooke's law was modified per the Lemaitre strain equivalence principle 31 , resulting in a new nonlinear statistical damage constitutive model for granite. The model's rationality and correctness were verified by comparing the theoretical values derived by the model with the experimental data.

Test method
Test samples. The rock samples utilized in this test came from the guide tunnel of Shuangjiangkou Hydropower Station's main and auxiliary workshops in China. The tunnel was dry and the value of the in situ pore pressure was 0, so pore pressure can be neglected. It took 15 days from the tunnel excavation to the extracted samples to be sent to the laboratory. After the irregular samples were sent to the laboratory, the samples were processed into cylindrical specimens with a diameter of 50 mm and a height of 100 mm in accordance with the method suggested by the International Society for Rock Mechanics 43 , as shown in Fig. 1a and Table 1. The specimens were dried in a constant temperature ventilated drying oven at 105 °C for 24 h and subsequently cooled in a sealed desiccator, followed by a series of tests as described below.
Test equipment. The tests were conducted using the rock mechanics test equipment (MTS815 Flex Test GT) of Sichuan University's College of Water Resources and Hydropower Engineering, as illustrated in Fig. 1b. The major technical parameters of the test equipment are as follows: the maximum axial force is 4600 kN; the maximum confining pressure is 140 MPa; the operating temperature is 20-200 °C; the axial and circumferential extensometer resolutions are ± 4 mm and − 2.5 to 8 mm, respectively; and the measurement accuracy is 0.5%. Test scheme. This study included five groups of conventional tests of granite, with designed confining pressures of 1, 5, 10, 30, and 40 MPa; these pressures were chosen based on the stress state of the rock mass at different depths from the underground powerhouse to the cave wall in Shuangjiangkou.
The experimental procedure during each test was as follows 43 : (1) an external load was supplied to the predetermined hydrostatic pressure at a rate of 0.05 MPa/s and allowed to stabilize; (2) an axial load was applied at a rate of 0.5 MPa/s while maintaining a constant confining pressure until the rock sample failed at the end of the test; and (3) the maximum axial force on the sample was measured, as well as the corresponding confining pressure.

Mechanical characteristics
Total stress-strain curve. The volumetric strain ε v of a cylindrical sample exposed to axial loading with a confining pressure can be calculated with the recorded axial strain ε 1 and lateral strain ε 3 .
The axial, lateral, and calculated volumetric strains vs. the applied deviatoric stress (σ 1 -σ 3 ) can be plotted to follow the route of the failure of a rock sample. Figure 2 shows five groups of conventional triaxial compression test curves for granite with various confining pressures. The increase in confining pressure has the following effects on rocks; that is, the peak deviatoric stress, which represents compressive strength, increases; the elastic limit increases; the deformation corresponding to the peak stress increases; and the rock properties also change, with plasticity and strain hardening becoming more apparent.
Strength and confining pressure. Experiments have shown that confining pressure may greatly improve the strength of most rocks. The axial peak strength of rock continues to rise as the confining pressure increases. Equation (1) and Fig. 3 show the relationship between confining pressure (σ 3 ) and peak deviatoric stress (σ 1 -σ 3 ) p . www.nature.com/scientificreports/

Energy evolution and deformation characteristics
Analysis of energy dissipation mechanism. As a geological structural material, rock is composed of one or more minerals. Due to the different physical properties, diagenetic conditions and processes of the compositions, the internal structure of rocks is very different, showing heterogeneity. There are numerous randomly distributed defects in rock materials at the microscale, mesoscale and macroscale levels, which cause the rocks to become damaged bodies. On the microscale, all deviations from the ideal lattice structure in the crystal configuration are called crystal defects. Dislocation is one of the most important kinds of defects whose slip movement causes crystal plastic deformation. The study of crystal plasticity, damage, strength, and other structural sensitivity problems relies heavily on dislocation theory 44,45 . On the mesoscale, the defects are mainly microvoids and microcracks whose nucleation and movement cause the damage evolution of materials. On the macroscale, the defect is a macroscopic crack and hole, which is common to our naked eye, and the appearance of macrodefects represents the beginning of material rupture. At a larger observation scale, defects are characterized by geological faults and large tectonic belts. The energy absorbed by rock as a result of an external load can be split into two categories. Releasable elastic strain energy is held in one component, whereas dissipated energy is dissipated in the other. The irreversible alteration of the internal structure of rock materials is always followed by an energy dissipation process. The energy dissipation mechanism of rock materials is very complex and can be roughly divided into four types. Two  www.nature.com/scientificreports/ main energy dissipation mechanisms are crystal plastic deformation and damage evolution. From the perspective of macroscopic mechanical behaviour, these mechanisms correspond to irreversible deformation and stiffness reduction; from the perspective of the microscopic structure of crystal materials, plastic deformation corresponds to crystal dislocation, and damage evolution corresponds to the nucleation and movement of microcracks 46,47 . Crystal plastic deformation and damage evolution are two types of closely related and interactive mechanisms. An inelastic deformation process generally has both types of mechanisms, and pure plasticity or pure damage is extremely rare. From a physical view, the two dissipative processes are essentially different; hence, plasticity and damage should be regarded as two independent energy dissipation mechanisms. During the loading process, rock materials contain a large number of microvoids and microcracks, the nucleation and movement of which are the dominant aspects of the energy dissipation mechanism, while crystal dislocation is a minor aspect. The friction of the fracture surface is also an important energy dissipation mechanism for rock materials, which is often the dominant energy dissipation mechanism under compression and shear conditions. Dilatancy caused by the roughness of the fracture surface is also an important aspect of energy dissipation, and the rock will show obvious dilatancy in the failure stage.
Energy calculation principle. If the conventional triaxial compression test system of rock is considered a closed system and the pressure system has no loss of external heat transfer during the process of loading rock samples, then the total energy generated by external work on rock samples can be calculated using Eq. (2) 16 .
where U, U d and U e represent the total absorbed energy density, dissipated energy density and released elastic strain energy density, respectively. The denser dotted area under the stress-strain curve reflects dissipated energy, while the sparser dotted area shows releasable elastic strain energy contained in rocks, as shown in Fig. 5. Equation (3) gives the elastic strain energy density based on elastic theory.
where ν and E are Poisson's ratio and the elastic modulus, respectively, which can be calculated from the linear elastic section of the stress-strain curve, as shown in Table 1.
The total absorbed energy density in conventional triaxial tests can be calculated by Eq. (7) 48,49 . www.nature.com/scientificreports/ where U 0 is the energy density due to hydrostatic loading, U 1 is the energy density due to axial compression by σ 1 after hydrostatic pressure, and U 3 is the energy density due to lateral compression by σ 3 after hydrostatic pressure. U 0 , U 1 , and U 3 can be determined by Eqs. (8)-(10) 48,49 .
where σ 1i , ε 1i , σ 3i , and ε 3i are the axial stress, axial strain, lateral stress, and lateral strain at point i on the stress-strain curve, respectively.
Energy analysis. Under external loading, the native microvoids and microcracks randomly distributed in rock first close; then a stress concentration will be produced at the tip. When the stress concentration exceeds the critical value, the crack of microcracks will occur and expand. With increasing load, microcracks will bifurcate and converge, forming a macroscopic fracture surface and even instability failure. During the process, energy is constantly dissipated. Figure 6 depicts the relationship between the deviatoric stress, energy density, AE event count, cumulative AE event count ratio and axial strain of Shuangjiangkou granite samples. Both the dissipation energy density and input energy density of Shuangjiangkou granite specimens increase as axial strain increases under various confining pressures. This indicates that the test system always inputs energy to the granite sample during the whole test process, and the energy is gradually dissipated due to plastic deformation (crystal dislocation), damage evolution (nucleation and movement of microcracks), friction of the fracture surface and dilatancy caused by the roughness of the fracture surface. When the elastic strain energy density curves of five groups of granite samples are compared to the stress-strain curves under the same confining pressure, it is discovered that the elastic strain energy density variation trend is very similar to that of deviatoric stress. The elastic strain energy density rises first and then falls, reaching its maximum value at peak stress with axial strain. This demonstrates that the prepeak period is mostly concerned with the storage of elastic strain energy, whereas the post-peak stage is primarily concerned with its release. At the elastic deformation stage, AE event rarely occurs. At the stage of steady crack development, AE event count increases. At the unstable crack development stage, AE event count increases sharply. When the stress reaches the peak, AE event intensifies. At the post-peak stage, AE event count reaches its peak. The variation trend of the cumulative AE event count ratio with axial strain is basically consistent with that of dissipated energy density, which indicates that the total energy absorbed by rock causes the AE event in the process of energy dissipation.
Energy dissipation ratio. The parameter energy dissipation ratio χ is introduced here, which can be calculated by Eq. (11). Figure 7 shows the variation trend of the energy dissipation ratio χ with axial strain under various confining pressures. The energy dissipation ratio has an obvious regular variation with axial strain, which first increases, then decreases and increases again. According to the five groups of curves (as shown in Fig. 7), curve fitting can be carried out for the five stages of energy dissipation ratio variation under different confining pressures; thus, the theoretical equation between the energy dissipation ratio and axial strain is obtained, as shown in Eq. (12) when the deviatoric stress (σ 1 -σ 3 ) is less than the crack initiation stress σ ci and in Eq. (13) when (σ 1 -σ 3 ) is equal or greater than σ ci . In order to have a simple functional relationship between energy dissipation ratio and axial strain, as well as between energy dissipation ratio and confining pressure, the coefficients are fitted with the confining pressure, as shown in Eqs. (14) and (15).
where r i and s i are the relevant parameters of the theoretical equation, which can be calculated by Eqs. (14) and (15), respectively.
where a i , b i , c i , d i , e i and f i are the relevant parameters of the theoretical equation, as shown in Table 2. www.nature.com/scientificreports/

Stress thresholds and deformation characteristics.
A large number of studies [11][12][13]50,51 show that a rock's stress-strain curve can be separated into five stages by using four threshold values (i.e., σ cc , σ ci , σ cd , and σ p ), as shown in Figs. 8 and 9. The ultimate stress (σ p ) is the peak point of the stress-strain curve, and the other three threshold values can be determined in many ways. Peng et al. 51 proposed a method to determine the crack    www.nature.com/scientificreports/ closure stress (σ cc ), which is easily programmed to make the determination more objective. Cai et al. 50 proposed a method to determine the damage stress (σ cd ) that has been proven to be relevant to jointed rock masses. A method called the lateral strain response (LSR) method used in this paper to determine the crack initiation stress (σ ci ) was introduced by Nicksiar and Martin 14 , which depends solely on the LSR and eliminates the need for user judgement.
A new method is proposed to determine σ cc and σ cd by using the relationship between the energy dissipation ratio and axial strain in this paper, as shown in Figs. 8 and 9. The stress at the curve's crest was determined to be σ cc , while the stress at the curve's trough was determined to be σ cd . Table 3 displays the stress thresholds of Shuangjiangkou granite samples subjected to various confining pressures.
As shown in Figs. 8 and 9, the deviatoric stress-strain curve of Shuangjiangkou granite obtained by conventional triaxial compression tests can be divided into five stages by four stress thresholds as well as the curve of the energy dissipation ratio and axial strain.
(1) Region I is the compaction stage of the rock's primary microcrack. The stress threshold corresponding to the end point is σ cc . The development process of rock deformation and failure is closely related to the changes in microcracks and microvoids that exist in the primary rocks. Under the condition of increasing stress, the primary microcracks and microvoids close under pressure at first, and the deviatoric stress-strain curve appears as an upper concave section. Most of the energy absorbed at the initial loading stage is stored in the rock, and a small part is dissipated due to rock compaction. However, the energy dissipation ratio increases rapidly, leading to an upwards trend of the curve. (2) Region II is the stage of linear elastic deformation. When the stress exceeds σ cc , the rock deformation enters the linear elastic stage. Although the axial strains and lateral strains also increase, microcracks and microvoids in the rock do not develop. Before the stress reaches σ ci , rock internal crystal dislocations are rare, and new microcrack nucleation and movement are also rare. Furthermore, no macroscopic surface fracture forms, so little energy is dissipated, and the majority of the energy absorbed is transformed into elastic strain energy. There is an obvious upper concave section in the elastic strain energy density curve. The energy dissipation ratio decreases, leading to a downwards trend of the curve. (3) Region III marks stable crack growth from the beginning of σ ci to the end of σ cd , which corresponds to long-term strength 15,52 . When the stress exceeds σ ci , microcracks and microvoids newly generated in the rock begin to develop, some microcracks connect to each other; the number of microcracks increases; some microcracks begin to merge, cross, and then gradually form irregular longitudinal macroscopic cracks; and the increasing rate of axial strain and lateral strain increases. The crystal dislocation and the nucleation and movement of microcracks in the rock develop stably. Crystal dislocation is the basis of plastic deformation whose macroscopic expression is irreversible deformation (Fig. 4). As the lattice array of microscopic structures within the rock slips, a row of side dislocations expands from the centre of the slip band. If the motion of these linear defects is intercepted by a granular boundary, then the number of dislocations will increase near the barrier, and the dislocation accumulation near the particles increases the local stress field so that the dislocations are forced to merge and form isolation zones. When a wide enough isolation zone forms, the adjacent planes are separated, and a microfracture nucleus forms. Once a microfracture occurs, the local stress field will change automatically. The increased stress causes further dislocation slip and new dislocation accumulation at the tips of microcracks. In other parts of the microcrack, stress will be released, allowing more dislocations to enter and widen the microcrack. Despite the fact that the energy dissipation ratio continues to fall, energy dissipation begins to rise, while elastic strain energy growth slows. Only a small portion of the absorbed energy has been dissipated, and the majority of it is still transforms into elastic strain energy. (4) Region IV is the stage of unstable crack growth. When the stress exceeds σ cd , the crystal dislocation emission and the nucleation and movement of microcracks in the rock begins to develop rapidly. Additionally, there are a greater number of microcracks, and the microcracks begin to merge, cross, gradually form irregular macroscopic crack. A macroscopic fracture surface begins to appear, and the rock enters the stage of rapid crack development, while work hardening can be observed. The speed of energy dissipation increases rapidly, and the curve of the energy dissipation ratio shows a clear upwards trend, while the growth rate of the elastic strain energy clearly slows. (5) Region V is the post-peak failure stage. As the stress exceeds the ultimate stress (σ p ), the rock failed rapidly.
Moreover, a number of relatively regular macroscopic cracks are formed inside the rock; these macroscopic cracks are connected to form one or more macroscopic fracture surfaces with a dip angle of less than 45° to the principal stress direction (Fig. 4). This behaviour is a typical failure under compressive stress and is generally described as compression-shear failure characterized by shear deformation along the fracture plane. Heat is generated by the intense friction of the macroscopic fracture surface; and the roughness of macroscopic fracture surfaces leads to obvious dilatancy. All of these factors cause the energy dissipation ratio to rapidly increase, and the curve shows a sharp upwards trend. After reaching the peak, the accumulated elastic strain energy is quickly released and then dissipated by converting it into various forms of energy, aggravating rock failure.

Nonlinear statistical damage constitutive model
Damage variable. Statistical approaches can be used to represent the properties of rock on the mesoscale.
Assume that a rock is made up of many meso-elements that are large enough to contain many microdefects. In comparison to the whole structure of the rock, the size of the meso-element is deemed to be appropriately minute. Therefore, the significant impact of a single microdefect can be ignored. This indicates that each meso- www.nature.com/scientificreports/ www.nature.com/scientificreports/ element can be considered a particle in rock in the context of continuum mechanics theory. It is then possible to describe rock damage and failure behaviour in terms of these meso-elements. Rock, which has obvious heterogeneity, is composed of one or several kinds of minerals. Native rock is filled with a random distribution of various microdefects, which are caused by nucleation and random movement of microcracks under loading, and the characteristics of the mechanical properties of microdefects show a random distribution, so the damage produced by the mesoscopic elements is also randomly distributed in the rock material. This indicates that rock's macroscopic characteristics are determined by their meso-element statistical properties. Equation (16) defines the damage variable that describes the granite damage characteristics.
where N D is the number of damaged meso-elements and N is the total number of meso-elements.
Considering that the damage of rock material in the loading process is a continuous process, the following simplified assumptions are made before modelling the rock: (1) rock can be regarded as an isotropic, homogeneous, continuous and brittle material in which there are microcracks and microvoids at the macroscale; (2) the meso-elements of rock obey Hooke's law before failure; that is, the meso-elements have linear elastic properties, and the elastic damage constitutive equation applies to each meso-element; and (3) the damage of rock is progressive, which is the gradual accumulation of the damage of the meso-elements.
As described in "Energy evolution and deformation characteristics", the energy dissipation mechanism of rock is very complex. Factors such as crystal dislocation, nucleation and movement of microcracks, frictional heat generation of macroscopic fracture surfaces, and dilatancy caused by roughness of macroscopic fracture surfaces all lead to the gradual dissipation of energy absorbed by rock from the outside, and energy dissipation occurs at the same time as rock damage evolution. The process of rock damage evolution causes energy dissipation. The study of rock damage evolution features from the standpoint of energy dissipation can, in essence, reflect the deformation and failure process of rock. The energy dissipation induced by damage is intrinsically random due to the random distribution of damage of meso-elements in rock material.
According to the previous analysis, the energy dissipation ratio has the characteristics of a random distribution. Suppose that the energy dissipation ratio (χ) follows a Weibull distribution 36 ; therefore, the probability density function of the energy dissipation ratio can be expressed by Eq. (17).
where P(χ), m, and k are functions of probability density, scale parameter and shape parameter, respectively.
The degree of energy dissipation of rock material when subjected to external force also indicates the degree of damage. Equation (18) can be used to express the total number of damaged meso-elements as the dissipated energy density grows.
Substituting Eq. (17)  Nonlinear statistical damage constitutive model. Strain and nominal stress can be directly measured experimentally in conventional triaxial compression tests of rock. According to the Lemaitre strain equivalence principle 31 , the damaged material's constitutive equation is the same as that of the original material, with the www.nature.com/scientificreports/ exception that the nominal stress is substituted with the effective stress. Equation (21) shows the relationship between the effective and nominal stresses.
where σ i and σ i are the nominal and effective stresses, respectively. The theory of deformation of porous materials containing a viscous fluid proposed by Biot 53 constituted an extension of Terzaghi's original one dimensional theory 54 and Hooke's law, which is shown in Eq. (22).
where p is pore pressure, 1/H is an additional physical constant which is a measure of the compressibility of the soild for a change in water pressure.
Substituting Eq. (21)  When the pore pressure is not considered, Y can be calculated by Eq. (30).
The parameters k and b can be calculated from a series of triaxial test data. The parameters related to the experiments without considering the influence of pore pressure are shown in Table 4.
Verification of the constitutive model. The comparison between the theoretical curves of the model proposed and the model from Cao 55 is shown in Fig. 10, in which also shown the experimental curves under various confining pressures. In the case of low confining pressure (as shown in Fig. 10a-c), the fitting effect of the model proposed in this paper is better than that of literature 55 at the pre-peak stage. In the case of high confining pressure (as shown in Fig. 10d,e), the fitting effect of the model proposed in literature 55 is slightly better at the elastic stage. However, the model proposed in this paper has better overall fitting effect. Due to the existence of confining pressure, the rock does not lose its ability to bear load immediately after reaching the peak strength, www.nature.com/scientificreports/ www.nature.com/scientificreports/ but there is a process of stress and strain readjustment, which leads to the unstable snap-back phenomenon of the test curve. Both models failed to simulate the phenomenon well, so more research is needed to solve the problem. The variation trend of the curves based on the theoretical model proposed is similar to that obtained from the experimental data, indicating that the theoretical constitutive model reflects the genuine constitutive relationship of granite, according to the analysis of the five sets of curves. The principles of elasticity, continuum mechanics, strain equivalence, and statistical analysis all serve as the theoretical foundation for the theoretical constitutive model. There are some theoretical assumptions and simplified processing when processing five groups of fixed confining pressure conventional triaxial test data through a series of mathematical analyses.

Conclusion
In this paper, MTS815 rock test equipment was used to carry out five groups of conventional triaxial tests on Shuangjiangkou granite, and the following conclusions were drawn after study: (1) Both the total input and dissipation energy density of Shuangjiangkou granite samples increase as the axial strain increases, whereas the trend of the elastic strain energy density is highly similar to that of the deviatoric stress. The variation of AE event count can predict the failure process of Shuangjiangkou granite. (2) The fitting formula of the energy dissipation ratio is given based on the curves of the energy dissipation ratio-strain under various confining pressures, and a new approach for calculating the crack closure stress and damage stress based on the energy dissipation ratio is provided. The deformation of brittle rock is inherently time dependent [56][57][58] , The model proposed didn't consider the effects of temperature and time, which requires more researches. The strength of granite given in this paper is short-term strength, however, the authors' another paper 59 studied the creep properties of granite collected from the same site. Therefore, the proposed model based on energy dissipation rate can provide reference for relevant scientific research, and when applying this model to guide engineering design, the influence of pore pressure must be considered, in addition, the effects of temperature and time should also be taken into account. www.nature.com/scientificreports/